from scipy.integrate import quad
result, error = quad(lambda x: x**2, -1, 1)
print(f"Quadrature result: {result}, Error: {error}")